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Abstract 

We perform a calculation of the relativistic corrections to the electromagnetic 
elastic form factors of the nucleon obtained with various Constituent Quark Mod- 
els. With respect to the non relativistic calculations a substantial improvement is 
obtained up to Q 2 ~ 2(GeV/c) 2 . 
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1 Introduction 



The non relativistic constituent quark models (CQM) have given good results in the study 
of the static properties of the nucleon [IJ ^J, like the baryon spectrum and the magnetic 
moments, and in a qualitative reproduction of the photocouplings p, [|, [5]. However, the 
standard CQ-Models are unable to reproduce the Q 2 behaviour of the electromagnetic 
form factors even in the low momentum transfer [0, ^, p|, |10j . 

The use of harmonic oscillator models give rise to form factors which decrease too fast 
with respect to the experimental data. Some improvement of this behaviour, specially 
in the case of the transition form factors, can be obtained by using more realistic wave 
functions |], |5|]. However, the problem of a reasonable description of the elastic and 
transition form factors of the nucleon in the framework of a Constituent Quark model is 
still open. 

The inclusion of relativistic effects is expected to be important in the description of the 
nucleon form factors. The structure of the electromagnetic current of a relativistic bound 
system still represents an unsolved problem. Much attention has been recently devoted 
to this problem, following substantially three main lines: the expansion of relativistic 
current operators in powers of the inverse quark mass, — , the evaluation of the current 
matrix elements in a light-cone approach and the expansion of the full relativistic current 
matrix elements, again in powers of — . 

The first approach takes into account the relativistic effects in the electromagnetic 
operators JTT], |HJ, the baryon states being the standard CQM ones. This type of rela- 
tivistic correction includes the qq-pair contribution to two-body currents coming from the 



one-gluon exchange [|13|, [14]]. The numerical results show that these effects are significant 



but not sufficient to explain the data. 

There are many interesting results obtained in a light-cone approach f| [15], |16| , such 
as the fact that the relativistic corrections to the transition form factors are important at 
low Q 2 || and that the root mean square radius of the proton is increased || This 



method is very useful since it allows to perform calculations starting from non relativistic 
wave functions. However there are still some problems in extracting form factors from 
the evaluated current matrix elements. 

In this work we follow the third method which consists of expanding the current 
matrix elements in powers 

of ^ 03, m and 

we propose a simplified approach useful for a 



preliminary calculation of the relativistic corrections to the elastic electromagnetic form 
factors of the nucleon obtained starting from different Constituent Quark Models. The 
use of Lorentz boosts for the quark spinors ensures that the relation between the dynamic 
variables of the initial and final states is relativistically correct. On the other hand, we 
assume that the quark internal motion is well described by the standard non relativistic 
wave function. The current matrix elements are constructed with a quark current operator 
containing only one-body terms and no quark form factors are introduced. We point out 
that the non relativistic expansion of the matrix elements of the present work, up to 
order m~ 2 , is coincident with that given by standard procedures M, [HJ, |IBJ introduced 



for the few-nucleon systems and no approximation is done with respect to the momentum 
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transfer Q 2 dependence. 

In Sec. 2, we describe the evaluation of the current matrix elements arriving at simple 
analytical expressions for the form factors. In Sec. 3, we discuss the results obtained with 
various 3q-wave functions and make a comparison with the experimental data. A brief 
conclusion is given in Sec. 4. 



2 The current matrix elements 

For the study of the transition process between the initial (I) and final (F) states, we have 
to calculate the current matrix element 

i=i 

where is the e.m. current of the i-th quark. We choose the Breit frame and so the 

— * — * 

total initial and final tetramomenta Pi = (Ei,Pi), Pp = (E F ,P F ) are related by 



Pi = ~ Pf = ~ \ , Ej = E F = y M 2 + ^ = E (2) 

where q is the virtual photon momentum, Q 2 = <f and M is the nucleon mass. We 
denote with p* (i = 1,2,3) the quark tetramomenta in the nucleon rest frame and we 
introduce the relative three-momenta 

Pp = ^(Pl -P2), PX = + P 2 ~ 2 P3 ) ( 3 ) 

— * 

which are conjugated to the standard Jacobi coordinates p and A. The 3-quark state is 
assumed to be 

= f[ BiUiip^^fx) , (4) 
i=i 

where the Bi (i = 1,2,3) are the usual Dirac boost operators that transform the 
quark spinors Ui(p*) from the nucleon rest frame to the Breit one. The boosted spinors, 
ipi = BiUi(p*), have the covariant normalization 

t/jii/ji = 1 . (5) 

In Eq. (4) 4>(p~p,p~\) is the standard non relativistic 3q-wave function, where for simplicity 
we have omitted the spin and isospin variables. The final state is written in a similar way. 
The current operator of the i-th quark, j^i), has the form 

3p{i) = i/TT li \ — . ( 6 ) 
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where m is the quark mass and e » (e [) is the initial (final) quark energy in the Breit 
frame. The normalization factors . f^, . f^-, have been introduced in order to obtain for 

the current matrix elements the correct expantion in powers of — ( i.e. coincident with 



m 

what is usually quoted in the literature) as shown in ref. p0| . The quark energies e j, e • 
are then expressed in terms of the corresponding quantities in the nucleon rest frame by 
means of standard Lorentz transformations. 

Finally, we add a factor 2E to the matrix element of Eq. (1) in order to take into 
account the normalization of the total matrix element. 

Because of the antisymmetry of the 3q-states, we can substitute Yh=i 3iJJ) with 
3 J M (3). The interacting quark absorbs the photon threemomentum in the Breit frame 
and therefore taking into account the Lorentz boost on the 3-quark, we can write the 
momentum conservation as follows: 




p'x = v\ - y , V P = P P (7) 

where the apices refer to the final momenta. The resulting expression for the current 
matrix element is complicated because of the presence of non local terms coming from 
the momentum dependence and the calculation can be performed numerically. However, 
in order to arrive at a preliminary calculation of the relativistic corrections to the e.m. 
current, we introduce some simplified assumptions. 

First, consistently with the use of a non relativistic model for the internal nucleon 
dynamics we approximate the quark energies in the nucleon rest frame as e* ~ m. 
Furthermore, we perform an expansion keeping contributions up to the first order in the 
relative quark momenta, but we treat exactly the dependence on the momentum transfer 
q. To this end, we introduce in the current matrix element of Eq. (1) the variable i?\ that 
is related to p\ and p\ ' in the following way 




PX = 7? A + n\lo— Q ■ ( 9 ) 

From the previous equations one also has 

n x = ^(p' x +p x ) . (10) 

We expand the current matrix element of equation (1) by keeping up to the linear terms 
in 7f\ and p p . In correspondence, the electric and magnetic form factors have zero- and 
first-order contributions. 

The results for the zero-th order charge and magnetic form factors can be given in a 
simple analytical form 

Gf(Q 2 ) = ^(ts)%G%(Q^) (11) 
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GfhQ 2 ) = ^sft^GZiQ^) , (12) 

where G 1 ^, and G 1 ^ are the electric and magnetic form factors as given by the non 
relativistic quark model. The quantities ts, tj and g a 

1 r MQ\ 
ts = -—[E V s - —i-} , (13) 



Mm 1 ' E 12 

Mm 



Em + 



M Q 2 
6E 



(14) 



with 



* = l + s (15) 

^ = K + ^Q 2 ] 1/2 , (is) 



36# 2 

M 2 ( 



m = [m 2 + -^Q 2 ] 1/2 , (17) 



as multiplicative factors. 

The first order contribution to the charge density matrix element is essentially of 
spin-orbit nature and a non relativistic expansion up to order m~ 2 gives the sum of the 
standard and the anomalous spin-orbit terms. The first order term will be omitted in our 
calculations, since it gives a numerically negligible contribution for nucleon states which, 
according to the models we use, are mainly in S-wave. The first order corrections to the 
magnetic form factors are of two types, spin-orbit like and convective. The first one can 
be disregarded for the same reason quoted above for the charge form factor, while the 
convective part gives in any case a small contribution. 

Therefore, within these approximations, the relativistic corrections introduce two kinds 
of modifications with respect to the non relativistic treatment: a multiplicative factor 
coming from the expansion of the quark spinors and the argument of the non relativistic 
form factors, i.e. the momentum transfer squared Q 2 , being replaced by Q 2 ^-. 

The current matrix elements must satisfy the current conservation equation 

%J$i = , (18) 

and it is satisfied in our case since the Constituent Quark Models we have used are based 
on local interactions. 



3 Results and comparison with experimental data 

The form factors of Eqs. (11) and (12) can be calculated using as input the nucleon form 
factors obtained in a non relativistic quark model. We present the results for different 
choices of the quark interaction, namely the h.o. [0] (Fig. 1), the three-body force hyper- 
central potential |H| (Fig. 2) and an exactly solvable potential based on a hypercoulomb 
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interaction || ( Fig. 2). All these models have been used for the description of the spec- 
trum |], H and of the photocouplings || |], [|, • The three-body force approach has 
allowed also a systematic analysis of the transition form factors for the excitation of the 



baryon resonances |T0[ . All of them contain also a spin dependent (hyperfme) interaction, 
which is essential for the description of the N — A splitting and for the excitation of quite 
a few resonances. For the elastic form factors, the configuration mixing coming from the 
hyperfme interaction does not produce strong effects, apart from the neutron charge form 
factor, and we shall omit it here. 

In the h.o. case, the choice of the h.o. parameter a is crucial. There are many different 



values of a quoted in the literature, according to the quantities to be fitted [f22j. We report 
in Fig. 1 the results obtained with 1) a = 0.229 GeV, which gives the correct r.m.s. 
radius of the proton |1|, |23| without the relativistic corrections, and 2) a = 0.410 GeV, 



which is necessary in order to reproduce the photoexcitation of the Di 3 and F 15 resonances 
|3|, f|] and corresponds to a confinement radius of the order of 0.5 fm. We note that the 
relativistic corrections increase the r.m.s. with respect to the non relativistic calculation 
231] . In fact from Eq. (11) one gets, for the h.o. proton charge form factor, 



where M is the nucleon mass and m = 4f. In order to get the correct radius one should 
use a = 0.285 GeV, which however is not too different from the choice 1). The results 
of Fig. 1 show that the relativistic corrections improve the h.o. form factors, but the Q 2 
behaviour is still different from the experimental data. 

In Fig. 2 we give the form factors obtained starting from the non relativistic calculation 



performed with the three-body force potential of ref. [21|. This potential has the form 
V(x) = — - + b con f x, where x is the hyperradius x = V P 2 + ^ 2 an d the values 
of the parameters are r = 4.59 and b con f = 1.61 jmT 2 . It should be noted that with 
this choice of the parameters and the inclusion of the standard hyperfme interaction, the 



three-body force allows to describe consistently the non-strange baryon spectrum pi]] , the 
photocouplings and the electromagnetic transition form factors [|Kj . 

In Fig. 2 we give also the results for the solvable model of ref. |§. It is based on the 
hypercoulomb potential Vh yc {x) = — - , with r = 6.39, to which a small confinement 
term is added. The advantage of this potential is that the results can be given in analytical 
form. For instance, the proton charge form factor is given by 

Ge'(Q 2 ) = , ■ 2 l Q2 7T - (20) 
i 1 + 24 t 2 m 2 J 2 

which, at variance with the h.o. case, for large Q 2 has a power-law behaviour. 

From the analysis of the results of Figs. 1 and 2, one sees that in general the inclusion 
of relativistic corrections improves significantly the non relativistic calculations. The 
improvement at low Q 2 is related to the correct non relativistic limit of the current matrix 
elements. The improvement at higher Q 2 depends on the relation between the initial and 
final state variables and allows to keep exactly the Q 2 dependence of the form factors. 
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It should be noted that the simultaneous reproduction of the spectrum and the pho- 
tocouplings requires a confinement radius of the order of 0.5 fm [[|, f|, |21|, |j and the 



relativistic increase is quite beneficial, but it is still not sufficient to get nearer to the 
data. In particular there still remain problems in the Q 2 behaviour in the low and medium 
range. Similar problems are encountered also in the transition form factors, both in the 
relativistic || and in the non relativistic [|K], |9| calculations, and so one can think that not 
only the relativistic corrections are responsible for the discrepancies between the CQM 
calculations and the experimental data. As already noted elsewhere |H| |^, ||], some fun- 
damental dynamical mechanism (effective at large distance, which means at low Q 2 )is still 
lacking, such as the explicit inclusion of quark-antiquark pairs both in the baryon states 
and in the electromagnetic transition operator. 

4 Conclusions 

We have calculated the relativistic corrections to the elastic nucleon form factors in a 
simplified and preliminary approach which leads to simple analytical expressions. We 
have used as input different Constituent Quark Models, namely the harmonic oscillator 



ly], the hypercentral model of ref. and the analytical model of ref. || showing that 
in all the three models the relativistic corrections are important, since they bring the non 
relativistic calculations nearer to the data, but still they are not sufficient. The persisting 
discrepancy may be an indication that further degrees of freedom (qq-pairs, gluons) should 
be included in the CQM in a more explicit way. 
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Figure 1: The charge (a) and magnetic (b) form factor of the proton and (c) the magnetic 
form factor of the neutron. The curves are the h.o. calculations using a = 0.229 GeV 
(dotted and dot-dashed) or a = 0.410 GeV (dashed and full). The dashed and dotted 
curves are the non relativistic calculations, the full and dot-dashed are the corresponding 
relativistic ones, obtained from Eqs. (11) and (12). The experimental data are taken 
from the compilation of ref. [15]. 
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Figure 2: The charge (a) and magnetic (b) form factor of the proton and (c) the magnetic 
form factor of the neutron. The curves obtained using the model of ref. [21] are the non 
relativistic (dashed) and relativistic calculations (full). The curves obtained using the 
model of ref. [9] are the non relativistic (dotted) and relativistic calculations (dot-dashed). 
The data are the same as in Fig. 1. 
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